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IS.  ABSTRACT 


This  report  presents  a  quantitative  measure  of  the  effects  of  energy  absorber  stroke  length  upon 
cumulative  probability  of  injury.  The  response  of  a  seated  subject  to  a  vertical  deceleration  pulse 
representing  a  helicopter  crash  was  calculated  for  a  given  seat  and  clothing  weight.  The  seat-man 
system  was  analytically  supported  by  a  "square  wave"  energy  absorber  selected  to  generate  no 
greater  than  23  G  for  a  5th  percentile  man.  The  calculated  response,  DRI,  provides  a  probability 
of  spinal  injury  for  the  seated  occupant.  Since  all  parameters  necessary  for  response  computations 
had  known  statistical  properties,  it  was  possible  to  calculate  the  joint  probability  of  injury  for  a 
particular  deceleration  pulse  and  subject  weight.  By  calculating  the  statistical  values  for  many 
deceleration  and  weight  combinations,  sufficient  to  represent  the  total  population  of  both,  a  cumu¬ 
lative  probability  of  injury  was  generated.  Stroke  length  of  the  energy  absorber  required  for  each 
combination  of  deceleration  and  weight  was  calculated.  By  examining  the  effects  of  a  limited 
stroke  length,  it  was  possible  to  generate  a  curve  of  stroke  length  available  versus  cumulative 
probability  of  injury.  The  curve  indicates  that  for  a  realistic  stroke  length  (12  inches),  the  cumu¬ 
lative  probability  of  injury  is  0.119.  By  doubling  the  stroke  available  or  by  halving  it,  the 
cumulative  probability  is  decreased  or  increased  by  7  percent.  Comparisons  with  previously  reported 
data  indicate  that  the  injury  potential  of  the  square  wave  is  significantly  higher  than  is  theoretically 
achievable. 
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This  report  was  prepared  by  Beta  Industries,  Inc. ,  under  the  terms 
of  Contract  DAAJ02-70-C-0055. 

The  purpose  of  the  effort  described  herein  was  to  determine  the  cu¬ 
mulative  probability  of  injury  to  a  helicopter  occupant  seated  in  an 
armored  crash-force-attenuating  seat  when  subjected  to  a  range  of 
possible  survivable  vertical  crash  forces.  The  scope  of  the  study 
was  limited  by  the  following:  (1)  only  one  energy  absorber  force- 
deflection  curve  was  used,  but  stroking  distance  was  varied;  (2)  the 
5th  through  95th  percentile  population  of  clothed  and  equipped  Army 
aviators  was  studied;  (3)  seat  bucket  weight  was  constant;  (4)  seat 
cushion  characteristics  were  not  varied;  and  (5)  only  vertical  crash 
forces  were  considered.  Crash  pulse  data  and  energy  absorber 
performance  criteria  were  extracted  from  USAAVLABS  Technical 
Report  70-22,  "Crash  Survival  Design  Guide." 

As  a  result  of  this  effort,  a  plot  of  cumulative  probability  of  injury 
versus  energy  absorber  limit  stroke  length  is  presented.  Correla¬ 
tion  with  previous  energy  absorber  efforts  indicates  that  for  the 
seat-man  system  studied,  probability  of  occupant  injury  may  be  re¬ 
duced  by  (1)  replacing  the  "square  wave"  energy  absorber  with  an 
energy  absorber  having  a  "notched"  force-deflection  curve,  (2) 
lowering  the  "square  wave"  energy  absorber  limit  load  below  that 
used  in  this  study,  or  (3)  using  an  energy  absorber  with  multiple 
load  limit  settings  as  discussed  in  USAAMRDL  Technical  Report 
71-22. 


The  results  of  this  effort  were  incorporated  in  the  recent  revision 
of  the  "Crash  Survival  Design  Guide,"  USAAMRDL  Technical 
Report  71-22,  revised  October  1971. 

The  Eustis  Directorate  technical  monitor  for  this  effort  was  Mr. 
G.  T.  Singley,  III,  of  the  Safety  and  Survivability  Division. 
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SUMMARY 


This  report  presents  a  quantitative  measure  of  the  effect  of  energy  - 
absorber  stroke  length  upon  cumulative  probability  of  injury.  If  a 
helicopter  crashes  with  an  energy  absorber  beneath  the  seated  occupant, 
the  acceleration  transmitted  to  him  will  generate  a  response  that  is 
indicative  of  the  probability  of  injury.  If  all  possible  crash  accelera¬ 
tion  environments  and  seat  occupant  weights  are  examined,  this  proba¬ 
bility  can  be  calculated  for  all  crashes.  Since  the  crash  environment 
can  be  described  statistically,  and  since  the  subject  weight  is  a  sta¬ 
tistical  parameter,  the  cumulative  probability  of  injury  ca.i  be  calcu¬ 
lated. 

A  digital  program  was  available  to  calculate  the  probability  of  injury 
for  a  subject  sitting  within  a  seat  supported  by  an  energy  absorber 
subjected  to  vertical  crash  impact  forces.  The  program  was  capable 
of  accepting  variable  weights  and  input  vertical  accelerations;  hence, 
it  was  possible  to  calculate  the  subject's  response  for  all  necessary 
combinations  of  inputs  and  weights.  The  data  were  collected  for  dis¬ 
crete  combinations  such  that  the  cumulative  probability  at  each  point 
could  be  calculated  as  the  product  of  the  individual  probabilities.  At 
each  point,  stroke  length  was  another  calculated  output.  By  summing 
the  probabilities  over  given  stroke  lengths,  it  is  then  possible  to  de¬ 
termine  the  effect  of  stroke  length  upon  cumulative  probability  of 
injury  for  a  given  seat  energy  absorber. 

The  curve  generated  indicates  that  for  a  realistic  stroke  length  (12 
inches),  the  cumulative  probability  of  injury  is  0.  119  By  doubling 
the  stroke  length,  the  cumulative  probability  is  decreased  by  7  per¬ 
cent  to  0.  111.  Halving  the  stroke  increases  the  probability  by  the 
same  value  to  0.  127.  Correlation  with  previously  developed  data 
contained  in  Reference  3  indicate  that  the  injury  potential  of  the 
square  wave  absorber  is  significantly  higher  than  is  theoretically 
achievable. 

The  data  presented  have  inherent  limitations  and  restrictions.  The 
results  were  computed  assuming  a  triangular  vertical  deceleration 
waveform  as  specified  in  the  Crash  Survival  Design  Guide.  ^  The 
subject  weights  were  based  upon  U.  S.  Army  Aviator  distribution, 
and  a  fixed  number  was  used  for  seat,  equipment,  and  clothing 
weight.  A  nonlinear  representation  of  a  polyurethane  foam  seat 
cushion  is  included  in  the  calculation.  The  DRI  concept  is  strictly 
applicable  to  spinal  injury  only  as  a  result  of  vertical  impacts. 
Additionally,  only  one  selected  energy  absorber  waveform  is  used 
for  all  cases  examined. 


FOREWORD 


The  statistical  study  reported  is  a  portion  of  an  overall  program  de¬ 
signed  to  develop  criteria  for  energy-abosrbing  landing  gear  as 
required  under  Contract  DAAJ02  -70-C-0055,  Project  IF  162203A529. 
This  particular  portion  was  concerned  with  the  effect  of  a  square -wave 
energy  absorption  waveform  upon  the  cumulative  probability  of  injury. 
The  computer  program  required  for  this  effort  was  available  and  the 
calculations  generated  data  indicative  of  human  response  to  input 
acceleration  waveform.  Hence,  it  was  possible  to  evaluate  a  crash 
condition  independent  of  development  of  a  particular  energy-absorbing 
landing  gear,  and  to  infer  seat-mounted  energy- absorbing  capability. 

Many  individuals  made  this  program  possible.  The  anthropometric 
data  were  provided  by  Mr.  George  Singley,  USAAMRDL;  the  DRI 
data  were  provided  by  Mr.  J.  W.  Brinkley,  U.  S.  Air  Force  Aerospace 
Medical  Research  Laboratories;  and  the  computer  time  was  supplied 
by  Mr.  R.  L.  Peterson,  Air  Force  Flight  Dynamics  Laboratory. 
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INTRODUCTION 


There  has  always  been  a  problem  in  determining  the  fox  ce  level  for  which  an 
energy  absorber  should  be  designed  in  order  to  protect  seated  subjects 
from  spinal  impact  injuries.  If  the  force  level  is  selected  assuming 
that  a  heavy  man  is  to  be  protected,  the  force  during  impact  ill  be  too 
great  for  a  lightweight  man.  If  the  force  level  is  selected  fox  the  small 
man,  then  the  large  man  will  require  a  greater  stroke  length.  If  the 
stroke  length  is  to  be  increased  to  avoid  bottoming,  for  which  acceler¬ 
ation  input  is  the  stroke  to  be  compatible?  The  inherent  difficulty  in 
determining  the  force  level  to  be  used  has  been  in  resolving  tie  amount 
of  injury  potential  the  device  possesses. 

A  recently  developed  computer  program  (3)  has  provided  the  necec  >ary 
tool  to  evaluate  the  injury  potential  of  a  particular  energy  absorber. 

The  analytical  model  is  shown  in  Figure  1.  The  program,  listed  in 
Appendix  III,  consists  of  nonlinear  ordinary  differential  equations  'ybich 
relate  the  forces  and  displacements  of  a  seated  subject,  structural 
seat,  and  energy  absorber.  Hence,  for  any  given  vertical  impulse  ii- 
put,  it  is  possible  to  calculate  the  response  of  the  seat  and  the  man. 

Since  there  exists  a  means  of  relating  seat  vertical  acceleration  to  in¬ 
jury,  an  additional  degree  of  freedom  is  added  to  the  "force”  model  of 
the  man  to  calculate  the  Dynamic  Response  Index  (DRI). 

The  desired  quantitative  measure  of  absorber  performance  is  obtained 
because  of  the  availability  of  statistical  data  relative  to  each  model 
parameter.  That  is,  the  crash  input  acceleration  has  a  statistical  dis¬ 
tribution  associated  with  it  to  relate  probability  of  occurrence  to 
deceleration  magnitude  (see  Reference  1).  Similarly,  the  crewmember 
weight  has  a  statistical  value  associated  with  it.  And  lastly,  a  given 
DRI  value  can  be  related  to  a  probability  of  injury  (5). 

With  both  deterministic  and  statistical  values  available,  bandwidth 
magnitudes  of  acceleration  and  subject  weights  can  be  used  as  input  to 
the  digital  routine  to  calculate  the  human  response.  Since  all  param¬ 
eters  are  statistically  independent,  it  is  possible  to  calculate  the  joint 
probability  of  injury  for  any  particular  crash  that  will  occur.  By 
summing  over  all  accelerations  and  weights,  it  is  possible  to  calculate 
one  number  which  indicates  the  cumulative  probability  of  injury  that 
exists  for  all  possible  vertical  accelerations  and  seated  subjects. 

Given  that  a  survivable  crash  occurs,  the  cumulative  probability  of 
injury  for  a  given  absorber  is  now  available. 

Another  parameter  automatically  calculated  is  the  stroke  length  re¬ 
quired  to  dissipate  the  input  energy.  If  a  stroke  length  is  selected  as 
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Figure  1.  Analytical  Representation  of  Seated  Subject. 
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being  the  practical  limit,  it  was  shown  (4)  that  the  impact  at  bottoming 
produces  a  DRI  increase  sufficient  to  cause  the  probability  of  injury  to 
approximate  unity.  Therefore,  once  the  array  of  responses  has  been 
generated,  the  variation  of  cumulative  probability  of  injury  is  evaluated 
by  setting  the  DRI  probability  to  unity  and  resumming  all  cell  values. 
This  was  accomplished  to  provide  a  quantitative  measure  of  the  effects 
of  changing  available  stroke  length,  for  a  given  energy  absorber. 
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CALCULATION  OF  INJURY  POTENT LAL 


The  purpose  of  this  research  was  to  calculate  the  variation  of  cumula¬ 
tive  probability  of  injury  for  changes  in  stroke  length  by  implementing 
an  existing  digital  program.  The  initial  steps  required  to  do  this  were: 

1.  Collect  the  statistical  data  required 

2.  Generate  the  computer  input  data 

3.  Determine  the  force  level  of  the  energy  absorber 

4.  Calculate  the  response  of  the  human  to  all  vertical  acceleration 

inputs  and  man-weight  variations 

5.  Generate  the  statistical  tables 

6.  Calculate  the  cumulative  probability  of  injury  as  a  function  of 

stroke  length 

These  are  depicted  in  Figure  2. 


STATISTICAL  DATA 


Statistical  data  were  collected  for  the  crash  deceleration  profile,  the 
seated  subject  weight,  and  the  probability  of  injury.  The  statistical 
data  for  the  crash  deceleration  pulses  were  obtained  from  the  "Crash 
Survival  Design  Guide",  USAAVLABS  TR  70-22  (1),  Figure  1-11 
(Distribution  of  Vertical  Impact  Accelerations  for  Rotary  -  and  Light 
Fixed-Wing  Aircraft).  The  distribution  of  personnel  weights  was  de¬ 
termined  from  "U.  S.  Army  Aviators,  1970"  (6)  Natick  Laboratories. 
Data  for  the  Dynamic  Response  Index  (DRI)  were  obtained  from  the 
Impact  Branch  of  the  Aerospace  Medical  Research  Laboratory,  Wright - 
Patterson  Air  Force  Base.  These  distributions  are  presented  as 
Figures  6,  7,  and  8  in  the  appendixes. 

COMPUTER  INPUT  DATA 


The  computer  program  (see  Appendix  II  )  used  sol  °s  for  the  injury 
potential  of  a  particular  seated  subject  responding  io  a  particular 
input.  That  is,  it  is  necessary  to  calculate  the  response  for  particular 
waveforms  and  particular  subject  and  seat  weights.  The  cumulative 
probability  of  injury  should  theoretically  be  calculated  from  continuous 
distributions  of  the  variable,  but  this  cannot  be  practically  achieved. 
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Hence,  it  is  necessary  to  utilize  discrete  values  of  the  variable  and 
establish  the  probability  that  exists  for  a  particular  bandwidth.  This 
was  done  for  both  the  deceleration  profile  and  the  subject  weight. 

The  helicopter  crash  deceleration  G  distribu  tion  (1)  used  was: 


Mean  G 

Probability 

Mean  G 

Probability 

1.80 

0. 4700 

18.  00 

0.  0133 

3.97 

0. 1700 

19.  87 

0.  0107 

6.  90 

0. 0900 

22.00 

0.  0080 

7.  93 

0.0700 

24.00 

0.  0080 

10.00 

0. 0400 

26.00 

0.  0080 

12.00 

0. 0400 

28.00 

0.  0080 

13.  75 

0. 0267 

30.  00 

0.  0080 

16.00 

0. 0133 

32.  00 

0.  0080 

34.  00 

0.  0080 

1. 0000 


The  weights 

of  men  (without  clothing  and  equipment 

weight)  used  were 

Weight 

Weight 

(lb) 

Probability 

(lb) 

Probability 

128.0 

0.  05 

172.  0 

0.05 

137.5 

0.  05 

175.  0 

0.05 

143.  5 

0.  05 

178.  5 

0.05 

148.5 

0.  05 

181.  5 

0.05 

152.5 

0.  05 

184.  5 

0.05 

156.0 

0.  05 

188.  0 

0.05 

159.  0 

0.  05 

193.  0 

0.  05 

162.  0 

0.  05 

198.  5 

0.05 

165.  5 

0.  05 

206.  5 

0.05 

169.0 

0.  05 

220.  0 

0.05 

The  equipment  and  clothing  weight  used  in  conjunction  with  the  above 
subject  weights  was  37.  4  pounds.  This  weight  included  all  clothing 
and  equipment  except  the  boots  and  was  comprised  of  such  items  as 
the  armored  vest  (15.  9  to  17.  7  pounds),  survival  vest  (5.  0  pounds), 
flight  helmet  (3.  3  pounds),  pistol  plus  twelve  rounds  (2.  5  pounds), 
and  clothing  (6.  0  pounds).  By  using  one  weight,  approximately  the 
median  based  upon  size,  the  variation  accepted  was  approximately 
2  pounds.  This  is  negligible  when  considering  the  other  quantitative 
aspects  that  are  within  the  overall  computation. 

Another  input  parameter  included  was  the  force -displacement  curve 
of  the  seat  cushion.  The  model  represents  a  1. 0  -inch  layer  of 
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polyurethane  foam  upon  a  0.25-inch  layer  of  closed  cell  ensolite  foam. 
The  force -deflection  curve  used  as  input  is  shown  in  Figure  9. 

The  helicopter  seat  bucket  weight  used  was  147  pounds.  This  was  pro¬ 
vided  as  being  indicative  of  the  weight  of  an  energy-absorbing  armored 
seat. 

ENERGY  ATTENUATION 


The  energy-absorber  force  level  was  dictated  by  the  procedures  out¬ 
lined  in  Reference  1,  Section  3.  3.  3.  That  is,  it  was  assumed  that  the 
attenuator  provided  a  square  wave  force-displacement  curve  and  that 
the  force  required  would  not  injure  the  5th  percentile  man.  Eighty 
percent  of  the  5th  percentile  man  with  clothing  and  equipment  is  136.  6 
pounds.  This  and  the  seat  weight  dictate  a  maximum  force  level  of 
6,  523  pounds  to  be  used  in  order  not  to  exceed  a  23  G  level.  Had  the 
heavier  man  been  chosen  with  a  29  G  limit,  the  force  level  would  have 
been  too  great.  The  force  deflection  curve  used  is  shown  below. 


FORCE 

(POUNDS) 


DISPLACEMENT 

(INCHES) 


CALCULATION  OF  INJURY  POTENTIAL 


The  analytical  model  used  in  calculating  the  response  of  a  seated  human 
to  spinal  impact  has  been  presented  and  documented  in  a  previous 
report  (3).  The  program  consists  of  a  set  of  nonlinear  ordinary 
differential  equations  that  are  solved  using  iterative  techniques.  For 
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any  given  input  vertical  acceleration  that  is  specified  as  a  function  of 
time,  the  time  dependent  forces,  accelerations,  and  displacements  are 
calculated.  The  human  response  is  calculated  in  two  ways  in  order  to 
evaluate  the  injury  level  (DR I).  The  equations  of  motion  are  written  for 
a  10  Hz  system  to  duplicate  man's  force  response.  This  yields  the 
correct  seat  acceleration  that  would  exist  if  a  human  were  subjected  to 
the  attenuated  acceleration  pulse.  The  seat  acceleration  is  then  used  as 
the  input  acceleration  to  the  DRI  model,  an  8  Hz  system.  The  output  of 
the  latter  model  is  then  the  desired  DRI  magnitude  which  is  directly 
related  to  probability  of  injury. 

For  each  of  the  discrete  values  of  input  acceleration  and  man  weight,  a 
calculation  was  made.  The  DRI  value  and  stroke  length  for  each  were 
recorded. 

STATISTICAL  TABLES 


Each  calculation  has  associated  with  it  three  probabilities.  The  pro¬ 
bability  of  being  within  that  weight  bandwidth  <pw>.  the  probability  of 

having  a  deceleration  within  the  given  bandwidth  (P^),  and  the  proba¬ 


bility  associated  with  the  calculated  DRI  (P  ).  The  computed  data 

DRI 

were  tabulated  in  the  format  shown  in  Figure  8.  The  upper  cell  of 
any  one  statistical  calculation  is  the  cumulative  probability  of  having 
the  acceleration  and  weight  magnitudes  occur  (P^)  (P^)  •  The 


middle  cell  is  the  probability  of  injury  (P^^),  and  the  bottom  cell  is 


the  cumulative  probability  of  all  three  variables 
P  =  <PA>  <PW>  (PDR1> 


The  upper  cell  of  the  second  table  is  the  magnitude  of  the  DRI,  and  the 
lower  is  the  stroke  length  required  to  absorb  all  of  the  input  energy. 
All  tables  developed  are  presented  in  Appendix  III. 


DEVELOPMENT  OF  CUMULATIVE  PROBABILITY  OF  INJURY 
VERSUS  STROKE  LENGTH 


The  generated  tables  consist  of  an  array  of  340  probability/stroke 
length  pairs.  If  the  probabilities  for  all  cells  are  summed,  the  cumu¬ 
lative  probability  is  that  associated  with  the  largest  stroke  length 
obtained.  The  heaviest  man  with  the  greatest  deceleration  requires 
nearly  28  inches  of  stroke.  If  one  absorber  of  this  length  were  pro¬ 
vided,  then  all  crash  occurrences  examined  would  be  protected  as 
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defined  by  Reference  1;  therefore,  the  cumulative  probability  of  injury 
(0.  109)  is  plotted  on  Figure  3  as  the  injury  potential  associated  with  a 
square  wave  absorber  of  at  least  28  inches'  stroke  length.  By  changing 
the  stroke  length,  it  is  apparent  that  the  injury  potential  changes.  If 
the  stroke  length  were  20  inches,  all  occurrences  examined  with  stroke 
lengths  calculated  at  greater  than  20  inches  would  be  instances  of 
"bottoming"  or  impact  of  seat  and  floor.  For  these  conditions  it  is 
reasonable  to  assume  that  the  "bottoming"  is  impulsive.  That  is,  the 
floor  is  stiff  enough  to  create  a  deceleration  pulse  such  that  the  body 
responds  as  a  function  of  the  velocity  change.  Since  the  DRI  is 
linearly  related  to  the  velocity  change  in  feet  per  second,  nearly  one 
to  one,  a  small  velocity  residue  of  1  foot  per  second  increases  the  DRI 
by  one  unit  (7).  A  one-unit  change  does  not  create  a  small  change  in  pro¬ 
bability  because  of  the  statistical  relation  of  DRI  to  probability.  Hence, 
even  for  a  small  velocity  residue  at  the  instant  of  "bottoming",  the 
probability  of  injury  approaches  unity. 

The  cumulative  probability  for  a  given  stroke  was  calculated  by  selecting 
a  particular  stroke  length,  determining  all  cells  where  that  stroke  was 
exceeded,  setting  the  DRI  probability  to  unity  for  those  cells,  generat¬ 
ing  a  new  joint  probability  table,  and  resumming  all  values.  Only  the 
DRI  probability  is  set  to  unity.  A  given  cell  has  a  probability  asso¬ 
ciated  with  the  weight  and  deceleration  magnitude  of  that  cell.  Hence, 
the  new  cell  generated  has  the  joint  probability  of  weight -deceleration 
as  its  final  value. 

The  process  discussed  above  was  used  rather  than  recalculating  on  the 
computer  the  DRI  magnitude  that  would  occur  for  a  restricted  stroke 
length.  This  was  used  before  and  verified  as  a  reasonable  approxima¬ 
tion  technique  in  previous  works  (3). 
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RESULTS 


The  data  of  Figure  3  indicate  that  there  would  be  a  possible  improve¬ 
ment  in  cumulative  injury  probability  of  23-1/2  percent  if  a  27.  7- inch 
stroke  were  possible.  For  a  more  realistic  stroke  length  of  12  inches, 
the  improvement  is  approximately  17  percent  over  the  no -energy - 
absorber  condition.  The  variation  of  percent  improvement  as  a  function 
of  stroke  length  is  shown  in  Figure  4. 

The  curve  shown  is  relatively  flat  in  that  although  the  maximum  stroke 
required  has  been  determined,  only  a  2  3-1 /2-percent  improvement  is 
possible.  The  reason  for  this  is  apparent  in  the  tabulated  data.  Al¬ 
though  the  energy  absorber  is  designed  not  to  exceed  23  G  for  the  5th 
percentile  man,  and  to  generate  only  approximately  18  G  for  the  heavy 
man,  the  response  of  the  human  is  much  greater  than  desired.  The 
95th  percentile  man  and  the  95th  percentile  crash  pulse  are  attenuated 
over  27.  7  inches.  The  response  of  the  man  is  a  DRI  of  20.  81,  or 
approximately  a  20 -percent  probability  of  injury  due  to  an  18  G  force 
level.  Similarly,  the  5th  percentile  man  attenuates  the  95th  percentile 
crash  in  17.  4  inches  while  reaching  a  DRI  of  28.  58.  This  is  the  result 
of  an  energy-absorber  designed  to  generate  a  23  G  force  on  a  rigid 
mass  subject. 

The  data  are  indicative  of  the  fact  that  the  human  is  not  a  "rigid  mass" 
for  the  environments  considered.  Admittedly,  the  seat  cushion  is  in¬ 
cluded  in  the  analysis.  However,  the  cushion  is  a  device  that  is  soft 
only  up  to  full  body  weight.  The  one  used  is  thin  and  is  essentially 
"bottomed  out"  prior  to  any  great  acceleration  input.  Hence,  the  data 
are  more  indicative  of  the  human  response  and  not  the  effects  of  the 
cushion. 

To  further  demonstrate  the  effects  of  body  dynamics,  the  curve  de¬ 
veloped  has  been  superimposed  on  data  developed  under  another 
program  (4).  The  curves  shown  in  Figure  5  are  for  "optimized", 
notched  energy  absorbers. 

The  "optimized"  force  levels  are  those  that  will  not  exceed  a  given 
DRI  for  the  mean  weight  subject,  a  115 -pound  seat  bucket,  and  the  same 
deceleration  pulse  used  for  the  square  wave  analysis.  The  "optimiza¬ 
tion"  process  is  based  upon  taking  advantage  of  the  dynamic  nature  of 
man's  force  and  injury  response.  This  creates  a  "notched"  energy- 
absorber  force -displacement  curve.  The  curves  are  for  a  theoretical 
waveform  and  as  such  as  idealistic  or  only  conceptually  attainable. 
However,  they  are  indicative  of  possible  goals.  One  particular  wave¬ 
form  is  shown  on  the  following  page. 
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Stroke  Length,  In. 


Figure  3.  Cumulative  Probability  of  Injury  Variation 
With  Energy -Absorber  Stroke  Length. 


11 


PERCENT 


2 


Figure  5.  Cumulative  Probability  of  Injury  Variation 
With  Comparative  Data. 
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DISPLACEMENT 

(INCHES) 


The  "optimized"  energy  absorber  curves  in  Figure  3  indicate  that  for 
a  12 -inch  stroke  it  is  possible  to  have  a  50 -percent  improvement  in 
cumulative  probability  of  injury,  compared  with  the  square  wave  ab¬ 
sorber,  by  designing  for  an  18  DRI  limit  for  the  mean  ( 170-pound)  man. 
Even  if  the  force  level  is  selected  to  develop  a  23  DRI,  instead  of  23  G, 
the  improvement  in  probability  from  square  wave  to  optimum  absorber 
would  be  30  percent. 
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CONCLUSIONS 


Several  conclusions  can  be  drawn  concerning  the  curve  developed 
(Figure  3): 

1.  The  cumulative  probability  of  spinal  injury  for  a  helicopter 

crash,  with  vertical  deceleration  only,  with  the  selected 
seat  cushion  and  a  12 -inch  stroke  "square  wave"  absorber 
defined  on  page  7,  is  0.  119.  Given  that  a  vertical  decelera¬ 
tion  survivable  helicopter  crash  occurs,  there  is  an  11.9 
percent  probability  of  spinal  injury.  This  is  approximately 
a  17  percent  improvement  over  the  no  energy  absorber  con¬ 
dition  (P  =  0.  124). 

2.  The  stroke  required  for  the  most  severe  condition  for  the 

seat/energy  absorber  system  assumed  is  27.  7  inches. 

3.  The  cumulative  probability  difference  between  a  no-energy- 

absorber  condition  and  the  maximum  stroke  condition 
(27.  7  in.  )  to  handle  the  heaviest  man  and  greatest  deceler¬ 
ation  is  23.  8  percent. 

Considering  the  relative  magnitudes  of  the  generated  curve  and  pre¬ 
viously  reported  data,  it  is  concluded  that: 

1.  A  cumulative  probability  of  injury  lower  than  that  created  by 

the  vertical  load  limiting  criteria  of  Reference  1  is  theo¬ 
retically  achievable. 

2.  For  a  12 -inch  stroke  length,  a  50 -percent,  improvement  in 

injury  probability  is  theoretically  possible  if  a  "notched" 
energy  absorber  is  used. 

3.  An  improvement  in  the  injury  probability  for  a  12 -inch  stroke 

length  "square  wave"  energy  absorber  is  possible  if  the  load 
level  setting  would  be  reduced  below  that  cited  in  1*. 

4.  An  improvement  in  the  injury  probability  for  a  12 -inch  stroke 

is  possible  if  an  energy  absorber  with  a  multiple  load  limit 
setting  capability  is  used  (2). 

*Note:  Reference  1  has  been  superseded  by  Reference  2,  which 
has  a  lower  load  limiter  setting  requirement. 
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RECOMMENDATIONS 


It  is  recommended  that: 

1.  Additional  effort  be  spent  in  deriving  an  acceptable  cumulative 

probability  of  injury. 

2.  The  energy -absorber  force  level  criteria  of  Reference  1  be 

revised  to  reflect  the  available  injury  potential  data.  Some 
level  of  DRI  should  be  selected  to  dictate  a  particular  design 
condition,  or  a  cumulative  probability  injury  level. 
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INPUT  DATA 
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CUMULATIVE  FREQUENCY  OF  OCCURRENCE  -  PERCENT 


Figure  6.  Distribution  of  Vertical  Impact  Accelerations 
for  Rotary-and  Light  Fixed-Wing  Aircraft. 
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Figure  7.  Weight  Distribution  of  U.  S.  Army  Aviators  (Ref.  3). 


19 


INJURY  PROBABILITY 


Figure  8.  Injury  Probability  Curve. 
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SEAT  DEFLECTION  -  INCHES 


Figure  9.  Seat  Cushion  Force-Displacement  Curve. 
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APPENDIX  II 


COMPUTER  PROGRAM  USERS  GUIDE 


The  computer  program  used  for  this  study  is  listed  on  the  following 
two  pages  and  is  discussed  in  the  subsequent  paragraphs  with  a 
complete  listing  of  the  terms  following  the  discussion. 

The  initial  cards  (CON)  establish  the  constants  to  be  used  in  specify  ¬ 
ing  the  structural  constants.  The  aircraft  floor  stiffness  is  (K6)  and 
damping  is  (C6).  The  impact  velocities  for  acceleration  waveforms 
are  Al,  A2,  and  A3. 

The  first  constant  function  cards,  CFN,  represent  the  acceleration 
waveform  used  as  an  input.  The  number  within  the  field  indicates 
the  number  of  time  data  points  within  the  field.  If  it  takes  10  points 
to  describe  the  acceleration  input,  CFN  (10.  0)  would  be  required. 

The  next  CFN  card  represents  the  energy  absorber  as  a  function  of 
displacement  and  a  series  of  three  points.  The  next  two  CFN  cards 
describe  the  seat  cushion  characteristics.  SC  is  the  force -deflection 
curve  for  the  seat  cushion  and  ISC  is  the  inverse  function  or  deflec¬ 
tion-force  curve.  This  curve  is  required  in  the  program  to  determine 
the  initial  deflection  of  the  seat  cushion. 

The  seat  weight  (W2  ),  aircraft  type  (AC),  energy  absorher  stroke  (SL) 
and  number  of  energy  absorbers  (NE)  are  read  in  as  constants.  If  it 
is  desired  to  vary  these  quantities,  they  can  easily  be  read  as  para¬ 
meters.  The  weight  of  the  seated  man  (WI)  and  the  input  reduction 
factor  vary  from  run  to  run  and  are  read  in  as  parameters. 

The  function  cards  establish  the  functional  relationship  between  the 
parameters  on  the  constant  function  cards.  That  is,  FUN  (AJ,  T)  in¬ 
dicates  that  the  data  on  the  CFN  card  AJ  is  a  function  of  time,  T.  The 
computer  will  establish  linear  relations  between  the  points  given  and 
interpolate  as  required. 

The  next  series  of  cards  merely  evaluate  the  stiffness  and  damping 
coefficients  for  the  models.  Since  the  man  models  are  fixed,  10  HZ 
and  ^  =0.  3  for  the  force  model  and  8  HZ  and  %  =0.  22  for  the  injury 
model;  the  coefficients  are  calculated  by  inserting  the  given  body 
weight. 

The  remainder  of  the  program  is  the  coding  of  the  differential  equa¬ 
tions  of  motion  previously  mentioned.  For  example,  the  differential 
equation  for  the  seat  is  initially  written: 

M2  *5!  =  Fea  +  Fst  +  W2  +  Fsc* 
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The  forces  in  this  equation  are  due  to  the  energy  absorber  (Fea),  struc¬ 
tural  loads  (Fg|).  inertial  force  (W2),  and  the  seat  cushion  (Fgc).  The 
seat  cushion  force  arises  from  the  relative  displacement  between  the 
man  and  the  seat  which  causes  the  seat  cushion  to  be  compressed. 

There  is  a  limit  to  the  amount  of  compression  the  seat  cushion  can 
withstand;and  when  it  is  reached, the  cushion  is  essentially  a  rigid  link 
between  the  man  and  the  seat  and  the  equation  of  motion  becomes 

OO  o  o 

M2X2  =  CiXj  +  KjXj  -  ClX2  -  KlX2  +  Fea  +  Fst  -  W2 

The  first  terms  are  directly  from  a  summation  for  forces.  The  last 
three  are  the  forces  due  to  the  energy  absorber,  forces  due  to  struc¬ 
tural  response,  and  the  inertial  term. 

As  the  energy  absorber  compresses,  the  force  generated  changes  with 
increasing  relative  displacement.  However,  it  cannot  displace  an  un¬ 
limited  amount.  If  an  energy  absorber  is  not  sufficient  to  absorb  the 
input  energy,  it  would  attempt  to  travel  a  distance  greater  than  prac¬ 
tically  available.  Therefore,  the  stroke  limit,  SL,  is  specified.  If 
the  device  attempts  to  exceed  this  limit.it  is  assumed  that  the  seat 
impacts  on  a  very  rigid  structure  having  stiffness  K6.  This  causes 
large-magnitude  and  high-frequency  "ringing''  to  occur,  and  structural 
damping  C6  is  included  to  limit  the  ringing. 

The  other  condition  is  that  of  seat  rebound.  If  the  absorber  dissipates 
the  energy  properly,  a  point  is  reached  where  the  relative  velocity 
between  the  floor  and  seat  goes  to  zero.  If  the  seat  and  man  want  to 
rebound,  the  energy  absorber  does  not  necessarily  respond  as  a  con¬ 
tinuous  element.  Since  rebound  data  is  not  available  for  all  types  of 
absorbers,  it  is  assumed  the  device  deforms  elastically  back  to  zero 
force  with  the  same  stiffness  as  the  original  elastic  deformation  of 
the  curve.  This  is  shown  below.  The  value  used  at  present  is  60,  000 
pounds  per  inch. 


FORCE 


DISPLACEMENT 


The  remainder  of  the  program  contains  the  logic  required  to  properly 
control  the  energy  absorber  "rebound"  and  the  structural  bottoming. 
The  structural  response  is  the  easiest  to  explain.  When  the  absorber 
deformation,  X6,  exceeds  the  selected  stroke,  SL,  the  displacement 
of  the  "rigid"  link  is  calculated  by  integrating  the  relative  velocity 
starting  at  the  time  when  X6-SL  is  greater  than  zero.  This  creates  a 
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force  STR  which  is  added  into  the  differential  equation.  When  the  de¬ 
formation  becomes  less  than  the  limiting  stroke,  the  integration  stops 
and  the  force  goes  to  zero.  The  seat  is  free  to  rebound. 

The  energy  absorber  action  is  more  complicated.  As  the  displace¬ 
ment,  X6,  increases;  the  relative  velocity,  IDX6,  eventually  slows  to 
zero.  At  this  point  the  absorber  wants  to  rebound,  but  not  along  the 
original  path.  Instead,  it  wants  to  unload  along  a  line  parallel  to  the 
initial  elastic  portion  of  the  force -displacement  curve.  There  the  last 
value  of  force  and  displacement  attained  at  zero  relative  velocity  must 
be  stored  as  the  initial  value  of  the  rebound  or  relaxation  portion  of  the 
curve.  This  is  the  purpose  of  the  "track  and  store"  operations  such 
as: 

DISP  =  TAS  (X6,  AN,  O.  ) 

FORED  =  TAS  (ENAC,  AN,  O.  ) 

From  these  initial  conditions  we  must  backtrack  along  the  stiffness 
curve  K7.  Therefore,  the  force  in  the  energy  absorber  is  the  stored 
value  minus  the  product  of  K7  and  X9.  When  the  absorber  reaches 
another  zero  relative  velocity  and  expands  as  far  as  it  will  go,  the  defor¬ 
mation  may  again  increase  along  the  elastic  curve  from  ihe  stored 
values: 


DIS2  =  TAS  (X6,  AB,  O.  ) 

FOR2  =  TAS  (ENAC,  AB,  O. ) 

The  occurrence  of  elastic  rebound  within  the  elastic  range,  even 
though  1DX6  goes  positive,  is  recognized  by  having  the  logical  test 
placed  upon  1DX7,  which  detects  the  first  change  in  relative  velocity 
as  positive  but  the  maximum  displacement  has  not  occurred  to  initiate 
1DX7.  After  the  peak  displacement  is  once  reached,  1DX7  has  some 
value  and  enables  the  computer  to  calculate  the  relaxation  effect 
rather  than  return  to  the  original  force  displacement  curve. 

The  output  of  the  program  can  be  selected  by  use  of  header  and  outs 
cards  in  the  formats  shown.  As  written,  the  program  prints  out  and 
plots  the  input  acceleration,  the  seat  acceleration,  the  energy  ab¬ 
sorber  crush  and  the  DRI.  These  provide  in  one  plot  the  input,  the 
seat  response,  energy  absorber  response,  and  the  human  physiological 
response. 
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COMPUTER  PROGRAM  GLOSSARY 


AC 

AF,  AH, 

DISP 

DRI 

EA 

ENAB 
ENAC 
FA,  HA, 

FOREA 

FORED 

FORS 

NE 

SL 

STR 


Aircraft  control  parameter.  This  is  input  to  the  pro¬ 
gram  to  select  either  the  fighter,  helicopter  or  transport 
acceleration  input  depending  upon  whether  or  not  the 
value  is  1,  2,  or  3  respectively. 

AT  Symbols  identifying  the  acceleration  data  points  fo* 
fighter,  helicopter,  and  transport. 

The  displacement  of  the  absorber  as  long  as  the  relative 
velocity  is  positive.  The  track  and  store  operation 
causes  DISP  to  track  X6  until  1DX6  goes  negative  and 
stores  that  value.  When  1DX6  goes  positive  again,  DISP 
tracks  until  another  velocity  reversal.  In  this  manner 
DISP  always  contains  the  initial  conditions  for  displace¬ 
ment  when  the  seat  attempts  to  compress  the  absorber. 

The  Dynamic  Response  Index  which  is  a  measure  of 
vertebral  fracture. 

Symbol  identifying  the  data  points  obtained  from  the 
energy  absorber  curve.  The  number  in  parentheses 
establishes  the  number  of  points  to  be  used  as  input. 

(Force  in  pounds  versus  time  in  seconds) 

Total  force  for  all  absorbers,  (pounds) 

The  force  developed  by  one  absorber  at  any  time,  (pounds) 

TA  Symbols  identifying  the  time  function  generated  by  the 
data  points  AF,  AH,  and  AT. 

Symbol  establishes  the  relation  between  the  energy  ab¬ 
sorber  curve  and  calculated  relative  displacement. 

The  force  the  absorber  develops  as  long  as  the  relative 
velocity  is  positive.  This  duplicates  the  actions  of  DISP. 

A  stored  force  at  the  beginning  of  expansion  or  com¬ 
pression. 

Number  of  energy  absorbers;  also  an  input. 

Stroke  limit  is  included  as  input,  (inches) 

The  force  developed  by  the  structure  at  impact  of  the 
seat  when  the  stroke  limit  is  exceeded. 
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Al,  A2,  A3  Velocities  of  the  fighter,  helicopter,  and  transport  air¬ 

craft  respectively.  These  are  used  as  the  initial  values 
of  the  integrators  of  the  program. 

Cl  The  damping  factor  of  the  10  HZ  man  model  (#/inch/sec. ) 

C3  The  damping  factor  of  the  8  HZ  injury  model  (# /inch/sec. ) 

C6  Damping  coefficient  of  structural  element  that  is  impacted 

when  seat  impacts.  (#/ inch/sec. ) 

DISC2  The  displacement  of  the  absorber  as  long  as  the  relative 

veloch"  *s  negative.  The  track  and  store  operation  is 
reversed  so  that  DIS2  will  always  contain  the  value  of 
absorber  displacement  as  the  absorber  attempts  to 
elongate.  This  provides  the  initial  conditions  for  absorber 
displacement  during  expansion. 

FOR2  The  value  of  the  force  in  the  energy  absorber  during  ex¬ 

pansion.  The  value  is  held  during  compression. 

K1  The  stiffness  of  the  10  HZ  man  model  (#/ inch) 

K3  The  stiffness  of  the  8  HZ  injury  model  (#/inch) 

K6  Stiffness  of  structural  element  that  is  impacted  when 

seat  impacts  after  exceeding  stroke  limit.  (#/inch) 

K7  Relaxation  stiffness  of  energy  absorber.  This  can  be 

an  assumed  equal  to  the  initial  stiffners  of  energy 
absorber  force-displacement  curve.  (#/inch) 

Ml  The  mass  of  the  10  HZ  man  model  (#/ inch/sec^) 

M2  The  mass  of  the  seat  (#/ inch/sec^) 

M3  The  mass  of  the  8  HZ  injury  model  (#/inch/sec^) 

X6  Relative  displacement  across  the  energy  absorber. 

2DX1,  1DX1,  The  acceleration,  velocity  and  displacement  of  the 
XI  force  model  mass.  Ml. 

2DX3,  1DX3,  The  acceleration,  velocity,  and  displacement  of  the 

X3  injury  model  mass,  M3. 

2DX4,  1DX4,  The  acceleration,  velocity,  and  displacement  of  the 

X4  crash  input. 
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Input  Reduction  Ratio. 


RA 

SC 

ISC 

FSC 

XSC 

1DX2, 

2DX51 

2DX52 

2DX5, 

X5 


Data  points  describing  the  force  deflection  curve  of 
the  seat  cushions. 

Inverse  of  SC  (deflection  force  data  points). 

Force  developed  by  the  seat  cushion  (pounds). 

Compression  of  seat  cushion  (inches). 

X2  The  velocity  and  displacement  of  the  seat  cushion  - 

man  interface. 

Acceleration  of  seat  after  the  seat  cushion  is  fully  com¬ 
pressed. 

Acceleration  of  seat  before  the  seat  cushion  is  fully 
compressed. 

1DX5,  The  acceleration,  velocity,  and  displacement  of  the 
seat. 
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APPENDIX  III 
COMPUTED  DATA 


Statistical  Tables  (Table  I) 


Acceleration 

Distribution 


P_A  x  P 
PDRI 

P 


Weight  Distribution 


W 


P^  is  the  probability  of  the  acceleration  occurrence. 
F*DRI  is  the  probability  of  injury  associated  with  the 
calculated  DRI . 

P  is  the  joint  probability  of  all  variables. 

P.rr  is  the  probability  of  the  weight  occurrence, 

w 


Acceleration 

Distribution 
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TABLE  I.  STATISTICAL  DATA 
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Table  I.  Continued 
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Table  I.  Continued 
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Table  I.  Continued 
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Table  I.  Continued 
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1'ABLE  II.  RESPONSE  DATA 
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Table  II.  Continued 
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Table  II.  Continued 
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Table  II.  Continued 
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Table  II.  Concluded 
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